Tracking arboviruses, their transmission vectors and potential hosts by nanopore sequencing of mosquitoes.

The risk to human health from mosquito-borne viruses such as dengue, chikungunya and yellow fever is increasing due to increased human expansion, deforestation and climate change. To anticipate and predict the spread and transmission of mosquito-borne viruses, a better understanding of the transmission cycle in mosquito populations is needed. We present a pathogen-agnostic combined sequencing protocol for identifying vectors, viral pathogens and their hosts or reservoirs using portable Oxford Nanopore sequencing. Using mosquitoes collected in São Paulo, Brazil, we extracted RNA for virus identification and DNA for blood meal and mosquito identification. Mosquitoes and blood meals were identified by comparing cytochrome c oxidase I (COI) sequences against a curated Barcode of Life Data System (BOLD). Viruses were identified using the SMART-9N protocol, which allows amplified DNA to be prepared with native barcoding for nanopore sequencing. Kraken 2 was employed to detect viral pathogens and Minimap2 and BOLD identified the contents of the blood meal. Due to the high similarity of some species, mosquito identification was conducted using blast after generation of consensus COI sequences using RACON polishing. This protocol can simultaneously uncover viral diversity, mosquito species and mosquito feeding habits. It also has the potential to increase understanding of mosquito genetic diversity and transmission dynamics of zoonotic mosquito-borne viruses.

Genomic epidemiology of human adenovirus F40 and F41 in coastal Kenya: A retrospective hospital-based surveillance study (2013-2022).

Human enteric adenovirus species F (HAdV-F) is a leading cause of childhood diarrhoeal deaths. The genomic analysis would be key to understanding transmission dynamics, potential drivers of disease severity, and vaccine development. However, currently, there are limited HAdV-F genomic data globally. Here, we sequenced and analysed HAdV-F from stool samples collected in coastal Kenya between 2013 and 2022. The samples were collected at Kilifi County Hospital in coastal Kenya from children <13 years of age who reported a history of three or more loose stools in the previous 24 hours. The genomes were analysed together with the data from the rest of the world by phylogenetic analysis and mutational profiling. Types and lineages were assigned based on phylogenetic clustering consistent with the previously described criteria and nomenclature. Participant clinical and demographic data were linked to genotypic data. Of ninety-one cases identified using real-time Polymerase Chain Reaction, eighty-eight near-complete genomes were assembled, and these were classified into HAdV-F40 (n = 41) and HAdV-F41 (n = 47). These types co-circulated throughout the study period. Three and four distinct lineages were observed for HAdV-F40 (Lineages 1-3) and HAdV-F41 (Lineages 1, 2A, 3A, 3C, and 3D). Types F40 and F41 coinfections were observed in five samples and F41 and B7 in one sample. Two children with F40 and 41 coinfections were also infected with rotavirus and had moderate and severe diseases as defined using the Vesikari Scoring System, respectively. Intratypic recombination was found in four HAdV-F40 sequences occurring between Lineages 1 and 3. None of the HAdV-F41 cases had jaundice. This study provides evidence of extensive genetic diversity, coinfections, and recombination within HAdV-F40 in a rural coastal Kenya that will inform public health policy, vaccine development that includes the locally circulating lineages, and molecular diagnostic assay development. We recommend future comprehensive studies elucidating on HAdV-F genetic diversity and immunity for rational vaccine development.

Coding-Complete Genome Sequences of Copenhagen and Copenhagen-Derived vP811 Strains of Vaccinia Virus Isolated from Cell Culture.

The coding-complete genomes of laboratory vaccinia virus strain Copenhagen and the Copenhagen-derived deletion strain, vP811, were determined by short-read sequencing. Relative to the NCBI reference genome M35027, seven common coding differences were revealed, including an intact copy of the vaccinia virus immunomodulator A46R in both Cop and vP811.

Multiplex MinION sequencing suggests enteric adenovirus F41 genetic diversity comparable to pre-COVID-19 era.

Human adenovirus F41 causes acute gastroenteritis in children, and has recently been associated with an apparent increase in paediatric hepatitis of unknown aetiology in the UK, with further cases reported in multiple countries. Relatively little is known about the genetic diversity of adenovirus F41 in UK children; and it is unclear what, if any, impact the COVID-19 pandemic has had on viral diversity in the UK. Methods that allow F41 to be sequenced from clinical samples without the need for viral culture are required to provide the genomic data to address these questions. Therefore, we evaluated an overlapping-amplicon method of sequencing adenovirus genomes from clinical samples using Oxford Nanopore technology. We applied this method to a small sample of adenovirus-species-F-positive extracts collected as part of standard care in the East of England region in January-May 2022. This method produced genomes with >75 % coverage in 13/22 samples and >50 % coverage in 19/22 samples. We identified two F41 lineages present in paediatric patients in the East of England in 2022. Where F41 genomes from paediatric hepatitis cases were available (n=2), these genomes fell within the diversity of F41 from the UK and continental Europe sequenced before and after the 2020-2021 phase of the COVID-19 pandemic. Our analyses suggest that overlapping amplicon sequencing is an appropriate method for generating F41 genomic data from high-virus-load clinical samples, and currently circulating F41 viral lineages were present in the UK and Europe before the COVID-19 pandemic.

Detection of human adenovirus F41 in wastewater and its relationship to clinical cases of acute hepatitis of unknown aetiology.

As of 8 July 2022, the World Health Organization (WHO) have reported 1010 probable cases of acute hepatitis of unknown aetiology in children worldwide, including approximately 250 cases in the United Kingdom (UK). Clinical presentations have often been severe, with liver transplantation a frequent clinical outcome. Human adenovirus F41 (HAdV-F41) has been detected in most children with acute hepatitis, but its role in the pathogenesis of this infection has yet to be established. Wastewater-based epidemiology (WBE) has become a well-established tool for monitoring the community spread of SARS-CoV-2, as well as other pathogens and chemicals. In this study, we adopted a WBE approach to monitoring levels of HAdV-F40/41 in wastewater before and during an acute hepatitis outbreak in Northern Ireland. We report increasing detection of HAdV-F40/41 in wastewater, concomitant with increasing numbers of clinical cases. Amplicon whole genome sequencing further classified the wastewater-derived HAdV as belonging to the F41 genotype which in turn was homologous to clinically derived sequences. We propose that WBE has the potential to inform community surveillance of HAdV-F41 and can further contribute to the ongoing global discussion supporting HAdV-F41 involvement in acute hepatitis cases.

Shotgun metagenomic sequencing of the first case of monkeypox virus in Brazil, 2022.

Monkeypox virus (MPXV), a zoonotic virus endemic to the African continent, has been reported in 33 non-endemic countries since May 2022. We report an almost complete genome of the first confirmed case of MPXV in Brazil. Shotgun metagenomic sequencing was completed in 18 hours, from DNA extraction to consensus sequence generation.

SARS-CoV-2 Testing in the Community: Testing Positive Samples with the TaqMan SARS-CoV-2 Mutation Panel To Find Variants in Real Time.

Genome sequencing is a powerful tool for identifying SARS-CoV-2 variant lineages; however, there can be limitations due to sequence dropout when used to identify specific key mutations. Recently, ThermoFisher Scientific has developed genotyping assays to help bridge the gap between testing capacity and sequencing capability to generate real-time genotyping results based on specific variants. Over a 6-week period during the months of April and May 2021, we set out to assess the ThermoFisher TaqMan mutation panel genotyping assay, initially for three mutations of concern and then for an additional two mutations of concern, against SARS-CoV-2-positive clinical samples and the corresponding COVID-19 Genomics UK Consortium (COG-UK) sequencing data. We demonstrate that genotyping is a powerful in-depth technique for identifying specific mutations, is an excellent complement to genome sequencing, and has real clinical health value potential, allowing laboratories to report and take action on variants of concern much more quickly.

Shotgun metagenomic sequencing of the first case of monkeypox virus in Brazil, 2022

ABSTRACT Monkeypox virus (MPXV), a zoonotic virus endemic to the African continent, has been reported in 33 non-endemic countries since May 2022. We report an almost complete genome of the first confirmed case of MPXV in Brazil. Shotgun metagenomic sequencing was completed in 18 hours, from DNA extraction to consensus sequence generation.

Resurgence of Ebola virus in 2021 in Guinea suggests a new paradigm for outbreaks.

Seven years after the declaration of the first epidemic of Ebola virus disease in Guinea, the country faced a new outbreak-between 14 February and 19 June 2021-near the epicentre of the previous epidemic1,2. Here we use next-generation sequencing to generate complete or near-complete genomes of Zaire ebolavirus from samples obtained from 12 different patients. These genomes form a well-supported phylogenetic cluster with genomes from the previous outbreak, which indicates that the new outbreak was not the result of a new spillover event from an animal reservoir. The 2021 lineage shows considerably lower divergence than would be expected during sustained human-to-human transmission, which suggests a persistent infection with reduced replication or a period of latency. The resurgence of Zaire ebolavirus from humans five years after the end of the previous outbreak of Ebola virus disease reinforces the need for long-term medical and social care for patients who survive the disease, to reduce the risk of re-emergence and to prevent further stigmatization.

A large epidemic of a necrotic skin infection in the Democratic Republic of São Tomé and Principe: an epidemiological study.

INTRODUCTION: In 2016-18, the Democratic Republic of São Tomé and Príncipe suffered a necrotic skin infection epidemic. METHODS: A surveillance system was established after increased hospitalisations for this infection. Microbiology results were available for samples analysed in December 2016 and March 2017 using whole genome sequencing and metagenomics. Negative binomial regression was used to study the association of weather conditions with monthly case counts in a time-series analysis. RESULTS: From October 2016 to October 2018, the epidemic cumulative attack rate was 1.5%. The first peak lasted 5 months, accounting for one-third of total cases. We could not conclusively identify the aetiological agent(s) due to the country's lack of microbiology capacity. Increased relative humidity was associated with increased monthly cases (incidence rate ratio (IRR) 1.05, 95% CI 1.02-1.09), and higher precipitation in the previous month with a higher number of cases in the following month (months with 0-49 mm rainfall compared with months with 50-149 mm and ≥150 mm: IRR 1.44, 95 % CI 1.13-1.78 and 1.50, 95% CI 1.12-1.99, respectively). DISCUSSION: This epidemic was favoured by increased relative humidity and precipitation, potentially contributing to community-based transmission of ubiquitous bacterial strains superinfecting skin wounds. FUNDING: World Health Organization Regional Office for Africa, Ministry of Health.

Rapid viral metagenomics using SMART-9N amplification and nanopore sequencing.

Emerging and re-emerging viruses are a global health concern. Genome sequencing as an approach for monitoring circulating viruses is currently hampered by complex and expensive methods. Untargeted, metagenomic nanopore sequencing can provide genomic information to identify pathogens, prepare for or even prevent outbreaks. SMART (Switching Mechanism at the 5' end of RNA Template) is a popular approach for RNA-Seq but most current methods rely on oligo-dT priming to target polyadenylated mRNA molecules. We have developed two random primed SMART-Seq approaches, a sequencing agnostic approach 'SMART-9N' and a version compatible rapid adapters  available from Oxford Nanopore Technologies 'Rapid SMART-9N'. The methods were developed using viral isolates, clinical samples, and compared to a gold-standard amplicon-based method. From a Zika virus isolate the SMART-9N approach recovered 10kb of the 10.8kb RNA genome in a single nanopore read. We also obtained full genome coverage at a high depth coverage using the Rapid SMART-9N, which takes only 10 minutes and costs up to 45% less than other methods. We found the limits of detection of these methods to be 6 focus forming units (FFU)/mL with 99.02% and 87.58% genome coverage for SMART-9N and Rapid SMART-9N respectively. Yellow fever virus plasma samples and SARS-CoV-2 nasopharyngeal samples previously confirmed by RT-qPCR with a broad range of Ct-values were selected for validation. Both methods produced greater genome coverage when compared to the multiplex PCR approach and we obtained the longest single read of this study (18.5 kb) with a SARS-CoV-2 clinical sample, 60% of the virus genome using the Rapid SMART-9N method. This work demonstrates that SMART-9N and Rapid SMART-9N are sensitive, low input, and long-read compatible alternatives for RNA virus detection and genome sequencing and Rapid SMART-9N improves the cost, time, and complexity of laboratory work.

Improvements to the ARTIC multiplex PCR method for SARS-CoV-2 genome sequencing using nanopore.

Genome sequencing has been widely deployed to study the evolution of SARS-CoV-2 with more than 90,000 genome sequences uploaded to the GISAID database. We published a method for SARS-CoV-2 genome sequencing (https://www.protocols.io/view/ncov-2019-sequencing-protocol-bbmuik6w) online on January 22, 2020. This approach has rapidly become the most popular method for sequencing SARS-CoV-2 due to its simplicity and cost-effectiveness. Here we present improvements to the original protocol: i) an updated primer scheme with 22 additional primers to improve genome coverage, ii) a streamlined library preparation workflow which improves demultiplexing rate for up to 96 samples and reduces hands-on time by several hours and iii) cost savings which bring the reagent cost down to £10 per sample making it practical for individual labs to sequence thousands of SARS-CoV-2 genomes to support national and international genomic epidemiology efforts.

Genomic Surveillance of Yellow Fever Virus Epizootic in São Paulo, Brazil, 2016 - 2018.

São Paulo, a densely inhabited state in southeast Brazil that contains the fourth most populated city in the world, recently experienced its largest yellow fever virus (YFV) outbreak in decades. YFV does not normally circulate extensively in São Paulo, so most people were unvaccinated when the outbreak began. Surveillance in non-human primates (NHPs) is important for determining the magnitude and geographic extent of an epizootic, thereby helping to evaluate the risk of YFV spillover to humans. Data from infected NHPs can give more accurate insights into YFV spread than when using data from human cases alone. To contextualise human cases, identify epizootic foci and uncover the rate and direction of YFV spread in São Paulo, we generated and analysed virus genomic data and epizootic case data from NHPs in São Paulo. We report the occurrence of three spatiotemporally distinct phases of the outbreak in São Paulo prior to February 2018. We generated 51 new virus genomes from YFV positive cases identified in 23 different municipalities in São Paulo, mostly sampled from NHPs between October 2016 and January 2018. Although we observe substantial heterogeneity in lineage dispersal velocities between phylogenetic branches, continuous phylogeographic analyses of generated YFV genomes suggest that YFV lineages spread in São Paulo at a mean rate of approximately 1km per day during all phases of the outbreak. Viral lineages from the first epizootic phase in northern São Paulo subsequently dispersed towards the south of the state to cause the second and third epizootic phases there. This alters our understanding of how YFV was introduced into the densely populated south of São Paulo state. Our results shed light on the sylvatic transmission of YFV in highly fragmented forested regions in São Paulo state and highlight the importance of continued surveillance of zoonotic pathogens in sentinel species.

Coast-to-coast spread of SARS-CoV-2 in the United States revealed by genomic epidemiology.

Since its emergence and detection in Wuhan, China in late 2019, the novel coronavirus SARS-CoV-2 has spread to nearly every country around the world, resulting in hundreds of thousands of infections to date. The virus was first detected in the Pacific Northwest region of the United States in January, 2020, with subsequent COVID-19 outbreaks detected in all 50 states by early March. To uncover the sources of SARS-CoV-2 introductions and patterns of spread within the U.S., we sequenced nine viral genomes from early reported COVID-19 patients in Connecticut. Our phylogenetic analysis places the majority of these genomes with viruses sequenced from Washington state. By coupling our genomic data with domestic and international travel patterns, we show that early SARS-CoV-2 transmission in Connecticut was likely driven by domestic introductions. Moreover, the risk of domestic importation to Connecticut exceeded that of international importation by mid-March regardless of our estimated impacts of federal travel restrictions. This study provides evidence for widespread, sustained transmission of SARS-CoV-2 within the U.S. and highlights the critical need for local surveillance.

Importation and early local transmission of COVID-19 in Brazil, 2020.

We conducted the genome sequencing and analysis of the first confirmed COVID-19 infections in Brazil. Rapid sequencing coupled with phylogenetic analyses in the context of travel history corroborate multiple independent importations from Italy and local spread during the initial stage of COVID-19 transmission in Brazil.

Coast-to-Coast Spread of SARS-CoV-2 during the Early Epidemic in the United States.

The novel coronavirus SARS-CoV-2 was first detected in the Pacific Northwest region of the United States in January 2020, with subsequent COVID-19 outbreaks detected in all 50 states by early March. To uncover the sources of SARS-CoV-2 introductions and patterns of spread within the United States, we sequenced nine viral genomes from early reported COVID-19 patients in Connecticut. Our phylogenetic analysis places the majority of these genomes with viruses sequenced from Washington state. By coupling our genomic data with domestic and international travel patterns, we show that early SARS-CoV-2 transmission in Connecticut was likely driven by domestic introductions. Moreover, the risk of domestic importation to Connecticut exceeded that of international importation by mid-March regardless of our estimated effects of federal travel restrictions. This study provides evidence of widespread sustained transmission of SARS-CoV-2 within the United States and highlights the critical need for local surveillance.

Genomic and Epidemiological Surveillance of Zika Virus in the Amazon Region.

Zika virus (ZIKV) has caused an explosive epidemic linked to severe clinical outcomes in the Americas. As of June 2018, 4,929 ZIKV suspected infections and 46 congenital syndrome cases had been reported in Manaus, Amazonas, Brazil. Although Manaus is a key demographic hub in the Amazon region, little is known about the ZIKV epidemic there, in terms of both transmission and viral genetic diversity. Using portable virus genome sequencing, we generated 59 ZIKV genomes in Manaus. Phylogenetic analyses indicated multiple introductions of ZIKV from northeastern Brazil to Manaus. Spatial genomic analysis of virus movement among six areas in Manaus suggested that populous northern neighborhoods acted as sources of virus transmission to other neighborhoods. Our study revealed how the ZIKV epidemic was ignited and maintained within the largest urban metropolis in the Amazon. These results might contribute to improving the public health response to outbreaks in Brazil.

Author Correction: Nanopore native RNA sequencing of a human poly(A) transcriptome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Importation and early local transmission of COVID-19 in Brazil, 2020

We conducted the genome sequencing and analysis of the first confirmed COVID-19 infections in Brazil. Rapid sequencing coupled with phylogenetic analyses in the context of travel history corroborate multiple independent importations from Italy and local spread during the initial stage of COVID-19 transmission in Brazil. (AU)

Nanopore native RNA sequencing of a human poly(A) transcriptome.

High-throughput complementary DNA sequencing technologies have advanced our understanding of transcriptome complexity and regulation. However, these methods lose information contained in biological RNA because the copied reads are often short and modifications are not retained. We address these limitations using a native poly(A) RNA sequencing strategy developed by Oxford Nanopore Technologies. Our study generated 9.9 million aligned sequence reads for the human cell line GM12878, using thirty MinION flow cells at six institutions. These native RNA reads had a median length of 771 bases, and a maximum aligned length of over 21,000 bases. Mitochondrial poly(A) reads provided an internal measure of read-length quality. We combined these long nanopore reads with higher accuracy short-reads and annotated GM12878 promoter regions to identify 33,984 plausible RNA isoforms. We describe strategies for assessing 3' poly(A) tail length, base modifications and transcript haplotypes.